Effective field theory for S — 1 quantum nematic 
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For S = 1 system with general isotropic nearest-neighbor exchange, we derive the low-energy 
description of the spin nematic phase in terms of the RP 2 nonlinear a-model. In one dimension, 
quantum fluctuations destroy long-range nematic (quadrupolar) ordering, leading to the formation 
of a gapped spin liquid state being an analog of the Haldane phase for a spin nematic. Nematic 
analog of the Belavin-Polyakov instanton with tt2 topological charge 1/2 is constructed. In two 
dimensions the long-range order is destroyed by thermal fluctuations and at finite temperature the 
system is in a renormalized classical regime. Behavior in external magnetic field is discussed. 
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Low-dimensional spin systems have been attracting 
permanent attention of researchers over more than half 
a century. A rich palette of their physical properties de- 
termined by the essential role played by quantum fluctu- 
ations makes them a very attractive playground for test- 
ing various theoretical concepts. In the last two decades, 
this interest has got a considerable impact, motivated 
particularly by the increasing availability of quasi-low- 
dimensional magnetic materials. A number of exotic 
"quantum spin liquid" states has been discovered, the 
most wide known example being the famous Haldane 
phase in integer-spin antiferromagnetic (AF) chainsi 

A generic example of the Haldane phase is the isotropic 
Heisenberg spin-1 AF chain. However, the most general 
isotropic exchange interaction for spin S = 1 includes 
biquadratic terms as well, which naturally leads to the 
model described by the following Hamiltonian: 

H= cos 9 (S n -S n+S )+ smO (S n -S n+S ) 2 , (1) 

<n<5> 

where S n are spin-1 operators at the lattice site n, 
and summation over the nearest neighbors is implied. 
There are indications*** that moderate biquadratic ex- 
change is present in the quasi-one-dimensional compound 
LiVGe2C>6. The points 8 — it and 8 — corre- 
spond to the Heisenberg ferro- and antiferromagnet, re- 
spectively. In one dimension (ID), the model Q is 
studied rather extensively, and a number of analytical 
and numerical results for several particular cases are 
available, 1 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 : 19 It is firmly es- 
tablished that the Haldane phase with a finite spectral 
gap occupies the interval —ir/4 < < 7r/4, and the fer- 
romagnetic state is stable for ir/2 < 8 < 5ir/4, while 
8 = 57r/4 is an SU(3) symmetric point with highly de- 
generate ground stated 

Exact solution is availabloi* for the Uimin-Lai-Suther- 
land (ULS) point 8 — ir/4 which has SU(3) symmetry. 
The ULS point was shown** to mark the Berezinskii- 
Kosterlitz-Thouless (BKT) transition from the massive 
Haldane phase into a massless phase occupying the inter- 



val 7r/4 < 8 < ir/2 between the Haldane and ferromag- 
netic phase; this is supported by numerical studies** 3 - 

The properties of the remaining region between the 
Haldane and ferromagnetic phase are more controver- 
sial. The other Haldane phase boundary 8 = —ir/4 
corresponds to the exactly solvable Takhtajan-Babujian 
modelji 5 - the transition at 6 = —ir/4 is of the Ising 
type and the ground state at 8 < —n/A is spon- 
taneously dimerized with a finite gap to the lowest 
excitationsi**** 7 * 8 *** 3 " The dimerized phase extends at 
least up to and over the point 8 = —ir/2 which has a 
twofold degenerate ground state and finite gapi 16 i 17 i 18 

ChubukoAiii used the Holstein-Primakoff-type bosonic 
representation of spin-1 operators 1 ' " based on the quadru- 
polar ordered "spin nematic" reference state with (S) — 



0, {Si. 



1, {S'l — 0), and suggested, on the basis of 



the renormalization group (RG) arguments, that the re- 
gion with 8 > 57r/4 is a disordered nematic phase. Early 
numerical studies^* seemed to have ruled out this possi- 
bility, and a common belief now—-*- is that the dimerized 
phase extends all the way up to the ferromagnetic phase, 
i.e., that it exists in the entire interval 57r/4 < 8 < 7ir/4. 
However, recent numerical results** indicate that the 
dimerized phase ends at certain 8 C > 57r/4, casting doubt 
on the conclusion reached nearly a decade ago. 

The aim of the present paper is to show that the low- 
energy dynamics of the model for 8 > 57r/4 can be 
effectively described by the nonlinear a model for a unit 
director field (i.e., a unit vector whose opposite directions 
are physically identical). The coupling constant becomes 
small in the vicinity of the ferromagnetic phase boundary 
8 = 5ir/4. This formulation allows one to establish many 
properties of the nematic phase by using extensive results 
available for the standard (vector field) 0(3) nonlinear a 
model. We also study the effect of external magnetic 
field, which can be easily incorporated in our formalism. 
We argue that in ID the ground state is disordered, in 
a complete analogy with the Haldane phase in case of 
the vector 0(3) model, and its elementary excitation is 
a massive triplet. In 2D, long-range nematic order exists 
only at T = 0. An explicit solution for the Belavin- 
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Polyakov instanton with half-integer charge in a (1+1)- 
dimcnsional isotropic nematic is presented. 

We start by introducing the following set of coherent 
states for S = 1: 

\u,v) = ^2(uj + ivj)\tj), je(x,y,z), (2) 
i 

where \tj) are three "cartesian" spin-1 states: 

\±) = T (1/V2)(\t x )±i\t y )), \0) = \t t ). (3) 

The coherent state is characterized by vectors u and v 
satisfying the normalization constraint u 2 + v 2 = 1. The 
freedom to choose an overall phase factor can be fixed by 
setting u ■ v = 0. It is easy to check that the resolution 
of identity / T>(u,v)\u,v)(u,v\ = 1 holds. 

In what follows we are interested in the region around 
6 = 5tt /A, hence it is convenient to use the notation 



B 



cos#=— Ji, sin# = — J 2 



J 2 > Ji > 



Using the states © , one can construct the coherent state 
path integral, and the effective Lagrangian will have the 
form 

L e $ = -2h v n ■ d t u n - (hn t n + s), (4a) 

n <nS> 

where the average of the local Hamiltonian for two neigh- 
boring sites 1 and 2 is, up to a constant, given by 

(^12) = -4Ji{(mi ■ u 2 )(y 1 ■ v 2 ) - («i • v 2 ){v 1 ■ u 2 )} 

- J 2 {(ui - u 2 -vi- v 2 ) 2 + (ui ■ v 2 + Vi ■ u 2 ) 2 } 

— B ■ (ui x vi + u 2 x v 2 ). (4b) 

Here we have included the Zeeman term — B ■ S n 
describing external magnetic field B. 

Assuming a uniform ground state and minimizing (H) , 
one arrives at the mean-field phase diagram shown in 
Fig. ^a): at zero field the ferromagnetic phase with 
u = v = l/v2 is stable for J 2 < J\, and at J 2 = J\ one 
has a degenerate first-order transition into the nematic 
phase with v = and parallel alignment of u (actually, u 
and v can be used on equal terms, and we just voluntarily 
choose v to be zero in the ground state). Note that vector 
11 is in this case a director since u and —u correspond to 
the same state. In external field the system acquires fi- 
nite magnetization (S) — B/[Z(J 2 — Ji)], where Z is the 
coordination number of the lattice, while nematic order 
persists in plane perpendicular to B. The magnetization 
increases with the field, and at B — Z(J 2 — J\) the ne- 
matic undergoes a second order phase transition into the 
phase with fully saturated magnetic moment. 

Our next aim is to study how the above classical mean- 
field picture changes due to quantum or thermal fluctu- 
ations. We pass to the continuum limit in Q, viewing u 
and v as smooth field distributions. From the mean-field 
solution one may assume that v <C u, and from the form 
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FIG. 1: Schematic T — phase diagram of the model (0 
in the vicinity of the critical point J 2 = Ji: (a) in dimension 
D > 2; (b) in one dimension. 



of the Lagrangian l|4a() it is clear that v plays the role of 
momentum conjugate to u, so that v will be eventually 
proportional to the time derivative of u (later this will 
be checked in a self-consistent way). We will thus keep 
only terms up to the second order in v, and derivatives of 
v will be neglected. Doing so, one obtains the following 
continuum version of the Lagrangian: 

L[u,v] = Vq 1 J d D r{ - 2hv ■ d t u- 2Z{J 2 - ,h)u 2 v 2 

z 

+ 2v(Bxu)-(J 2 /2)J2((S a - V)m) 2 }, (5) 

a=l 

where Vb is the volume of the elementary cell of the D- 
dimensional lattice, S a are vectors describing the position 
of Z nearest neighbors with respect to a given lattice site, 
and constraints u 2 +v 2 = 1, u v = are implied. In what 
follows, we will assume for simplicity that the lattice is 
hypercubic, then Z = 2D, Vq = a D , and (Sk ■ V) = aVk, 
k = 1 . . . D, where a is the lattice constant. 

The "slave" variable v under the assumption v <C u 
can be integrated out, yielding 



[2Z{J 2 - Ji)]" 1 {(J5 x u) - hd t u}. 



(6) 



Substituting 10 back into JSJ gives the following effective 
Lagrangian depending on u only: 



L c ff — 



d D 
a 



B x u\ 2 



D 
k=l 



(V fe «) 2 }, 
(7) 
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where c = [2ZJ 2 (J 2 — Ji)] a/% is the characteristic 
limiting velocity, and u now should be considered as a 
unit vector, u 2 = 1. Note that according to a change 
of sign of u automatically means a sign change for v, 
so that u is in this approximation completely equiva- 
lent to — u. The above description is valid at the en- 
ergy scales E < Eq = 2Z(J 2 — Ji). One readily ob- 
serves that is nothing but the Lagrangian of the well- 
known nonlinear a model used as the effective theory 
for antiferromagnetsiS^ (without the topological term). 
Even the additional terms in the second line of Eq. (7J , 
describing the effect of the external magnetic field, are 
identical to those appearing in the a model for antifer- 
romagnets. Thus, the low-energy dynamics of model 
in the nematic phase is similar to the dynamics of an 
antifcrromagnet, with the only yet important difference 
that instead of the unit vector of sublattice magnetization 
one now has the nematic director u: the order parameter 
space is RP 2 instead of S 2 . The cr-model formulation, in 
contrast to the spinwave approach of Chubukov, 11 allows 
one to study full nonlinear dynamics of the problem. 

At zero field, one can rewrite the Lagrangian JJJ in 
a standard notation using dimensionless space-time vari- 
ables x = (xq,x), x — r/a, xq — ict/a. The effective 
Euclidean action takes the following compact form: 



(a) 



-4b _ 1 f ( du 

h 2g J \ dx^ 



7-D+l 



X, 



with the coupling constant g is defined as 
g = {Z(J 2 - J0/2J 2 } 1/2 . 



(8) 



(9) 



Note that smallness of the coupling constant does not 
require a large-.? approximation, and is controlled solely 
by the closeness to the ferromagnetic phase boundary. 

In one dimension (D = 1 ) continuous symmetry can- 
not be broken, and the ground state of © is disor- 
dered with exponentially decaying correlations. The cor- 
relation length £ for the usual 0(3) (vector) a model 
can be obtained within Polyakov's RG approach^ as 
£o(3) ~ ae 2 ^/ g . In the RP 2 cr-model, however, there 
is a rescaling in flow equations because of the change in 
the measure: the physical field is not it, but the bilinear 
projector R = u T ® u. The action can be rewritten as 



T ' = Tg J ^ d » R ' d * R ) d ° +1: 



(10) 



where (A, B) — tr(A T B) denotes the scalar product. The 
/9-function in the leading order is the same as for the 0(3) 
modelfSl f3(T) = — ^T 2 , with the trivially rescaled cou- 
pling constant T — 2g. Thus, for the correlation length 
in the RP 2 cr-model one obtains 



RP 2 



ae 



ir/g — „pT V /j 2 /(J2-Ji) 



(11) 



in agreement with Chubukov's one-loop RG resulliAi for 
interacting spin waves. The elementary excitation is a 
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FIG. 2: Belavin-Polyakov-type soliton in a spin nematic: (a) 
the solution is given by tci (z) in the upper half-plane (respec- 
tive to the line connecting points a and b), and by W2(z) in 
the lower half-plane, with the appropriate branch cuts shown 
as thick solid lines and numbers denoting the phase of w; (b) 
schematic view of the solution, nematic director u is shown as 
an ellipsoid in a projection onto the figure plane; white spot 
indicates the end of ellipsoid that is under the figure plane. 



massive spin-1 triplet, and the gap A = hc/£ opens up 
exponentially slow as one moves away from the phase 
transition point J 2 = J\. 



A~2[J 2 (J 2 - Ji)] 



i/2 p -t\AVW>- j i) 



(12) 



The RP 2 and 0(3) er-models are also different with re- 
spect to their topological excitations. In the 0(3) model 
there is a localized solution with nonzero 7r 2 topological 
charge Q — J d 2 £,e^u ■ (d^u x d v u), known as the 
Belavin-Polyakov instanton (BPI)i2& The simplest BPI 
with Q = 1 is described by id = (z — a)/(z — b), where 
the complex variable w(u) = (ui+iu 2 ) / (I — U3) is defined 
as a function of the complex coordinate z = x\ + ixo, and 
generally any analytical function w{z) yields a solution. 28 
The BPI action Abpi = AirhQ/g does not depend on 
the parameters a, b which can be interpreted as coordi- 
nates of elementary entities, merons, constituting a BPI. 
It was speculatedfSSi 2 ^ that the correlation length £0(3) <x 
e A B pi/2h j g j-giatgd to the concentration of merons. 

In the RP 2 case the director nature of the field makes 
possible BPI- type defects with half-integer Q, whose ac- 
tion is exactly one half of that for their 0(3) a model 
counterparts: Indeed, consider a solution of the form 

w = wi{z) = yfEf ■ F° r the 0(3) cr-model such a solu- 
tion would be invalid, because it has a branch cut. In a 
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nematic, however, w(u) and w(—u) are physically identi- 
cal, and the above solution can be matched with another 

one, w = 102 (z) = \J fE^j so that W2(u) = tOi(— u) on 
some line. This is easily achieved by choosing the cuts 
as shown in Fig. [21 This solution has Q = ±4, and its 
action is just one half of that for the 0(3) BPI. Curi- 
ously, this fact correlates with the extra factor | in the 
correlation length exponent (|1 1|) . 

In the RP 2 model there is another type of topologi- 
cal defects, disclinations, characterized by a -k\ topolog- 
ical charge (vorticity) q. It is argued that their presence 
could produce the BKT transition in the isotropic casem^ 
However, one can see that such a transition would occur 
above some critical value of the coupling gsKT of the or- 
der of 1, and as long as g = (1 — J\j J2) 1 / 2 <C 1, one may 
expect that the disclinations will be bound in pairs and 
their effect can be neglected. 

Our approach easily allows one to incorporate the ef- 
fect of an external magnetic field. Weak magnetic field 
B < A acts on the spectrum only in a trivial way (the 
Zeeman shift), but at B = A the gap closes and the 
system enters the critical phase with algebraically de- 
caying correlations, which can be characterized as the 
Tomonaga-Luttinger liquids^! The resulting phase dia- 
gram for ID case is shown in Fig. ^b). 

For D = 2, at zero temperature the ground state 
should have long-range nematic order, in agreement with 



recent numerical results^ as long as g is small compared 
to some finite value g c of the order of 1 which marks 
the transition into a quantum disordered phase. This 
latter transition is expected to be the same as a finite- 
temperature transition in the three-dimensional classical 
Lebwohl-Lasher model, which is the first order suppos- 
edly due to the effect of disclination lines ^ At T = 
the phase diagram in presence of magnetic field has the 
mean- field form of Fig. QJa) . At T 7^ and B = ne- 
matic order is destroyed by thermal fluctuations, with 
the correlation length £ ~ ae 27rj2 ^ T (note extra factor ^ 
in the exponent, as compared to the standard resuh™&). 

In summary, we have shown that the low-energy dy- 
namics of the bilinear-biquadratic 5=1 system JQ| for 
~ 57r/4 can be effectively mapped onto the RP 2 non- 
linear a model. We have argued that in one dimension 
this model exhibits a disordered nematic state, support- 
ing early proposition of Chubukov 11 and recent numerical 
results" against the commonly adoptedSLSSiSi point of 
view. Using parallels with the extensively studied vec- 
tor version of the cr-model, one can easily extract nec- 
essary information on the properties of nematic phase. 
An instanton solution of the Bclavin-Polyakov type with 
half-integer topological charge is presented. 
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